Momentum space properties from coordinate space electron density 
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Electron density and electron momentum density, while independently tractable experimentally, 
bear no direct connection without going through the many-electron wave function. However, in- 
voking a variant of the constrained-search formulation of density functional theory, we develop a 
general scheme (valid for arbitrary external potentials) yielding decent momentum space properties, 
starting exclusively from the coordinate- space electron density. Numerical illustration of the scheme 
is provided for the closed-shell atomic systems He, Be and Ne and for Is 1 2s 1 singlet electronic ex- 
cited state for Helium by calculating the Compton profiles and the (p n ) expectation values derived 
from given coordinate space electron densities. 
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For a quantum mechanical N-electron system such as 
an atom, molecule or a solid, the coordinate space one- 
electron density n(r) is derived from the pertinent config- 
uration space many-electron (antisymmetric) wave func- 
tion 9(fi, r 2l ^3, rjv) through the marginal distribu- 
tion n(r) = Nj |*(f, r 2 ,r 3 , . . . ,r N )\ 2 d 3 r 2 d 3 r 3 . . .d 3 r N 
(electron spins may also be included). Analogously, 
the one-electron momentum density 7(75) is obtained 
from the N-electron momentum space wave func- 
tion <fr(p\,p 2 ,p 3 , • ■ • ,Pn) via the reduction j(p) = 
N f\<I>(p,p 2 ,p 3 ,...,p N )\ 2 d 3 p 2 d 3 p 3 ...d 3 p N . The wave 
functions \& and $ in the complementary spaces are con- 
nected by a many-particle Fourier-Dirac transformation : 

$(pi,p 2 ,...,p N ) = 1 , / *(ri,f 2 ,...,fjv) 



(2tt) 3JV / 2 

EN _ _ 
i=i P3r3 d 3 r 1 d 3 r 2 ...d 3 r N . 



x e 
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(Hartree atomic units, viz. \e\ = l,m e = \,h = 1 have 
been used throughout). Experimentally, the coordinate- 
space density is tractable through a coherent and elastic 

x-ray scattering process |lj where the scattered intensity 

-> 2 

is directly proportional to |/(fc)| , f(k) being the form 
factor which is the Fourier transform of n{r). On the 
other hand, the electron momentum density ^{p) man- 
ifests itself more directly in terms of the (directional) 
Compton profile J(q) [2j obtained in an inelastic high 
energy (X-ray or 7-ray) Compton scattering process : 



dp x 



dPyl(Px,Py,q)- 
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For atomic and molecular systems in gaseous state, a 
spherically symmetric Compton profile results from the 
corresponding spherically averaged electron momentum 
density ~f sp h(p) = 37 / l{p)dVl p , whence 



Js P h(q) 



q>0 



= 2tt 



Js P h(p) P dp, 
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which is a monotonic decreasing function of q. Theoret- 
ically, the expressions in Eqs. J3J and J2J are essentially 
the impulse approximation (IA) profiles 3]. Equation 
© readily leads to an inverse relation "/(p) — > 
where here and henceforth, the subscript "sph" will be 
understood. 

It must be emphasized here that the mappings n — > ^ 
and 7 — > $ are both in general one-many and although 
for the ground state the former is unique 0, the ex- 
plicit prescription for the map is unknown. Hence while 
there exists a Fourier connection between $ < — ► <£, no 
such direct relation is possible between the densities n{f) 
and "f(p) in the two complementary spaces. On the ba- 
sis of quasi-classical phase-space considerations (akin to 
the Thomas-Fermi theory) jthere exists a procedure due 
to Burkhardt 0, Konya @, and Coulson and March 
0, called the BKCM scheme, to estimate electron mo- 
mentum density, given its position-space counterpart . 
However, this method is marred by its artifacts of a di- 
vergent 7(0) and a finite cut-off for . Incidentally, the 
so termed "Wigner function" 0, 0, [lij cannot represent 
a joint probability in phase space for not being strictly 
positive semidefinite. In their phase-space approach to 
density functional theory (DFT) Q , Parr et al. ^3 P re_ 
scribed a phase-space entropy maximization, imposing a 
given coordinate density and a given kinetic energy den- 
sity (at each point r) as constraints. This enabled them 
to obtain a positive semidefinite phase-space distribution 
through which momentum density could be extracted. 

Let us, however, pose a question : Given exclusively 
the electron density in coordinate space as a starting 
point {and no other information) , could one estimate the 
quantum-mechanical electron momentum density (and 
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hence the Compton profile)? 

It is the spirit underlying this letter to demonstrate 
an affirmative answer to the above question, within the 
density functional theory pertinent features of which have 
been as highlighted below. 

In their exciting work, Zhao and Parr (ZP) ^jj devel- 
oped a novel method to obtain the Kohn-Sham orbitals 
for a given co-ordinate space density. Their method is 
based on Levy's constrained search approach 0] that 
generalizes the Hohenberg - Kohn formulation of DFT. 
Constrained search approach which obtains the Kohn- 
Sham "wave function" , a single Slater determinant 
formed out of the lowest occupied orbitals of a local po- 
tential, by minimizing the non-interacting kinetic energy 



T s [n] = min < y D \f\^ D >, 



(4) 



where ty D =4> n(r), the given density; f is the N-electron 
kinetic energy operator and ty D is an antisymmetric, 
normalized (hence L 2 ) N-electron wave function of in- 
dependent electrons. ZP accomplished the search on the 
right side of Eq. Q through variation of the orbitals 
{tpi{r)}fLi yielding a density n(f). This density would 
equal the given density no(r) at every point r, via the 
minimization of the positive semidefinite functional 



c = f Kr) - n (r)][n(r') - n (r')} 



d 3 rd 3 r', (5) 

2\f-r'\ 

whose minimum value zero would be reached iff n(r) = 
rig (r) Vr . The minimization T s [n] + AC with respect to 
the orbitals {ipi} gives a set of Kohn-Sham like equa- 
tions : 



+ Xv{r) - Z/r ) iPi(r) = e^i(f), (i = 1, . . . , N) 

(6) 



where 



the sum 



\r—r' | 



\v(r) — Z/r 
' in the limit 



A 



Xv(r) 
oo,C - 
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yet the product AC remaining finite, gives the effective 
Kohn-Sham potential Rj. 

Morrison and Zhao 15j ap plie d this method to atoms 
while Ingamells and Handy extended the work to 
molecular systems. Recently, Harbola |J3 observed that 
the Zhao-Parr procedure could also be applied to obtain 
Kohn-Sham orbitals for an excited-state density, thereby 
demonstrating the general applicability of the method to 
ground- as well as excited-states. 

Our scheme to effect n(r) — > j(p) "transformation" 
is : Start from a given density no(r) => Obtain the 
Kohn - Sham orbitals ipi{r) via the Zhao - Parr pre- 
scription I ^iif) I = n{r) => Fourier transform 

— > 4>%ip) Obtain = J2°=i I <f>i(p) | 2 ; bence 
the Compton profile J(q) and other momentum expecta- 
tion values. This procedure thus starts from only a given 
rio(r) and estimates J(q) and (p n ) values. We illustrate 
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FIG. 1: Electron Momentum Density obtained within the 
present scheme 



this "r-density to momentum-space properties" (RDMP) 
scheme for the following atomic systems : 

Accurate coordinate space densities employed as start- 
ing points for ground state are due to Koga et al. 0] 
for He, Esquivel and Bunge 0] and Bunge and Esquivel 
[20| for Be and Ne respectively while the Coolidge and 
James [21| density was used for He excited state. The 
value of A (cf. eq. El ) was set to ^5000, leading to a 
sufficiently accurate self-consistent Kohn-Sham potential 
converged to five places. 

A collage of the jip) derived from n(f) under the 
present scheme is depicted in Fig. ^ While all the 7-plots 
compare extremely well with their HF-counterparts (not 
shown), the features of nonmonotonicity of j(p) for Neon 
is also reproduced (cf. Ref. [22J; Fig. 1), in conformity 
with the fact that atoms with their ultimate p-shells dou- 
bly occupied or more evince such nonmonotonicity . 

The Hartree-Fock data (zeroth-order correlated) for 
the wave-function are given for comparison and as a da- 
tum. Particularly striking is the fact that for excited He, 
the electron momentum density (EMD) is overwhelm- 
ingly preponderant around low momentum values along 
with a very sharp asymptotic fall-off. 

Tables [Q through HTH give the CP's via the RDMP 
scheme with their other theoretical and experimental 
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TABLE I: Compton Profiles for helium 



TABLE II: Compton Profiles for Ground State of Beryllium 
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comparators. Each CP is normalized as J °° J(q)dq = 
N/2, as is customary. There is a remarkable agreement 
between the CP's derived from coordinate space atomic 
electron densities within the impulse approximation and 
the Hartree - Fock (HF) or experimental CP's. From Ta- 
ble [J the considerable redistribution of electron momen- 
tum density of He excited state in comparison with its 
ground state counterpart is evident. Excitation brings 
in a slower decay in the coordinate space which is, by 
Fourier reciprocity, mapped on to a peak Compton pro- 
file value. As the transition Is 1 2s 1 — > Is 2 is dipolc- 
forbidden, this singlet He-excited state is long-lived (life 
time ~ y-jsec); while there have been no Compton pro- 
files (CP's) reported for this system, the present scheme 
accomplishes this. 

For Be and Ne (Tables |H] and 1111(1 ground states, the 
CP's via the present schemes give better overall agree- 
ment with their accurate correlated counterparts than do 
the HF-CP's. This is indicative of the correlation piece 
picked up by the Zhao-Parr scheme. 

To gauge the overall quality of the electron momen- 
tum density "derived" from the coordinate space den- 
sity, we have computed the (p n ) values under the RDMP 
scheme and have compared with those derived from ac- 
curate, correlated atomic wave functions 1281 . as well as 
from the near Hartree- Fock wave functions 23] . Table llVl 
illustrates that the present RDMP scheme successfully 
obtains the (p n ) values (for n = -2, -1, 1, 2, 3 and 4) 
in agreement with both their HF and correlated counter 
parts. For the case of He excited state, preponderance of 
electron momentum density around very small as well as 
large p values is conspicuous. 

Since the Kohn - Sham theory, in its very spirit, pro- 
vides an effective local one- body potential in which 
the mutually noninteracting electrons are immersed, 
the quantal exchange-correlation effects of the electron- 
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TABLE III: Compton Profiles Ground State of Neon 
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electron interactions (after filtering out the "classical" 
part J n(r')d 3 r' / \ f — r' | in v^A(r*)) embody a ki- 
netic energy like piece and a potential energy like piece. 
The kinetic piece arises out of the difference between the 
Junctionals T[n] and T s [n], whose values albeit known, 
their forms remain unknown. Lam and Platzman [30| 
and Tong and Lam |2o| imported the functional forms 
for a homogeneous electron gas and estimated the cor- 
rection to J(q) within the local density approximation. 
In the present case however, the difference as a func- 
tional cannot be isolated from the sum total Kohn-Sham 
exchange-correlation energy. 

To conclude, it is gratifying that the present method 
offers a general prescription to estimate quantal mo- 
mentum space properties starting from coordinate space 
n(r) alone, with no reference to the many- electron wave 
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TABLE IV: Various expectation values (p n ) for n = -2, -1, 1, 2, 3 and 4 
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function. This scheme could also be extended to solids 
wherein directional C omp ton profiles could be derived 
from an experimental |3l| three-dimensional co-ordinate 
space density. Note that this procedure is not limited 
by the form of the external binding potential endowing 
the present scheme with generality for any bound state 
problem. 
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